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Деякі проблеми обчислення трансформацій гравітаційного поля 


Точні трансформації гравітаційного поля в регіональному масштабі забезпечує розв'язання 
нелінійної граничної задачі Алексідзе для потенціалу сили тяжіння. Її заміна на відповідну задачу 
для рівняння Лапласа генерує похибки, оцінки яких наведені в статті. Чисельне моделювання 
останньої задачі виявило, що точність її розв'язання істотно залежить від міри обумовленості задачі, 
наближень напряму шуканого градієнта, параметрів чисельного методу та напряму зовнішньої нормалі, 
уздовж якої обчислюють похідні. 
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Методи трансформацій. Результатом гравіметричних зйомок є набір значень модуля градіє- 
нта потенціалу ий сили тяжіння (МГПСТ) |; - ягадиі. -и(5) на деякій замкнутій денній пове- 
рхні 5 , яка цілком охоплює тяжіючі маси в області С пошуку трансформант. Задача його пря- 
мих трансформацій нетривіальна через невідомий напрям вектора я. і не належить до ліній- 


них граничних задач математичної фізики. При наближеному розв'язанні цієї задачі її граничні 
дані можна задавати на незамкнутій, але досить великій поверхчхні (11. 
Загалом задачу трансформації значень гравітаційного поля слід вирішувати шляхом розв'я- 


зання граничної задачі безпосередньо для диференціального оператора А(а) сили тяжіння в 


області задання 8. Взагалі нелінійне рівняння для сили тяжіння має вигляд |2| 


А(з) - Дя(х, ) - а"(я,х, (ж, ) є м(х; ) х, ЄС, нзб8. (1) 
3 оре д 2, 
де а"(з,х,)- У сова, Асова, -у Рана -о геометричний фактор поля тяжіння, 
й іч іс МОХ; 5 
ух, ) 2 -Алу|вгай сх, ) - задана функція, пов'язана з розподілом мас, Є, - кут між напрямом 


вектора я|; та і-ї координатної осі. Це рівняння враховує напрям сили тяжіння і складне в пра- 
ктичній реалізації. Якщо вектор яЇ; у кожній точці простору спрямований в одну і ту ж точку 


Р(хо)), і - 13 (уздовж нормалі до поверхні еліпсоїда), рівняння (1) спрощується (ЗІ: 


Ї 





Один із практичних способів трансформації і перерахунку значень сили тяжіння - розв'язан- 
ня граничної задачі безпосередньо для значень МГПСТ |4), інший - розв'язання нелінійної зов- 
нішньої граничної задачі Алексідзе для потенціалу сили тяжіння: 





Лінеаризує цю задачу припущення щодо напрямку МГПСТ на поверхні 5, але воно вносить 
значні похибки, а визначення напряму - не менш громіздка задача, ніж визначення МГПСТ. Зо- 
крема, для зовнішньої тяжіючої точки отримано рівняння |3| 


З о ВЕНВИ 
зАсо8а, 3,8, з - Де со5 И, - 0, і - 1,3; 
ТИЙ ХХ; 


Ї 


які визначають модуль сили тяжіння Є за відомим розподілом кутів й, в просторі С або коси- 


нуси кутів 2, за відомими значеннями Є у просторі С (виникає невизначеність СО5, на о). 


Питання розв'язності і єдиності задачі (3) принципове для гравіметрії: воно відповідає, чи 
можливе високоточне однозначне визначення сили тяжіння в зовнішньому просторі за високо- 


точними абсолютними значеннями сили тяжіння на замкнутій поверхні 5, яка цілком охоплює 


усі маси, і які потрібні мінімальні додаткові дані. За одночасного визначення на поверхні 5 аб- 
солютних значень і напряму сили тяжіння її обчислення в зовнішньому просторі нескладне, але 
такі роботи дорогі і в масових розвідувальних роботах не проводяться. 

Розв'язують задачу (3) шляхом послідовних наближень. У плоскому випадку задача визна- 
чення гармонічної функції, модуль градієнта якої заданий на замкнутому контурі, не має єди- 
ного розв'язку (31. 


Якщо граничні значення у (5) і функція й залежить від двох змінних, гранична задача (3) 
стає плоскою 





Якщо 12) з и(х, чн ім(Х, ), і- 12 - довільна голоморфна в С функція, реальна частина 
якої співпадає з розв'язком граничної задачі (2), її уявна частина теж є розв'язком задачі (4). 


Співвідношення Коші-Рімана КК. --Ї, де ка ди Ре іа бу | ду кутові коефіцієн- 
дх, | ОХ, Ох ОХ, 


ти прямих, спрямованих уздовж градієнтів функцій и і у, є умовою їх перпендикулярності. 








Нехай у скінченному И - околі напряму Ї градієнта РР о одного з розв'язків гранич- 
ОЇ 


ної задачі (4) не існує іншого розв'язку задачі (4), тобто, не існує двох розв'язків задачі (4), 
напрями градієнтів яких складали б менший за // кут. Кут між градієнтами двох розв'язків зада- 
чі (4) - га АТ а розв'язки, що різняться між собою константою - тотожні. 


Для перевірки цього припущення чисельно змодельовано методом розкладення за неортого- 
нальними функціями |5) граничну задачу (4) за допомогою збіжного ітераційного процесу 


Ди" (х.)-0, Уа) уз); да (5) 
іі Х; 


І 


де и" - К-те наближення розв'язку граничної задачі (4), а коефіцієнти напрямних косинусів 


о, (шо )- о У и 


Похибки трансформацій. Похибка є заміни рівняння (2) на рівняння Лапласа Де (х,)- 0 на 





2 








отримані з нормування граничних умов задачі (4). 


сфері С радіусу а з початком координат у її центрі, в точці Х, - 0, і - Ї,3З дорівнює 


а" у, (б«є ча" у /6 (6) 


де за умови, що область С дотична до поверхні Землі, нормальне гравітаційне поле: 


. 77 2.981000 мГл 
у тах віх //г є Зтах (о) піт аз 6360? км 
| . 2-(978000--300а) мГл 
(6380- а) км! 


звідки, згідно (6), одержуємо при а - 100 км - 76 «є «81 мГал. При перерахунку сили тяжіння у 
вільному просторі за формулою Пуассона похибка може сягати 100 мГал на великих площах 


У, я тій 28 (х, у «2тіп є (х, У пах | 


х:єС х|єС хієЄС 


земної поверхні. У плоскому випадку при х, - 0, 2-13 дає 51 «є 554 мГал. За (6), Є - ? -9 «0, 


тому дійсне значення сили тяжіння завжди менше, ніж одержане у припущенні її гармонічності. 
Порівняння похибок вказує, що рівняння сили тяжіння краще апроксимувати рівнянням (2), ніж 
Лапласа. 

У практиці гравіметрії при перерахунку в зовнішній простір розв'язують не внутрішню гранич- 
ну задачу, а зовнішню. Похибка перерахунку у зовнішньому просторі (заміну зовнішньої гранич- 


ної задачі Де (х.)- 22 (х.)/ г -0 на зовнішню задачу Да (х.)-0 за однакових граничних умов 





5 


реа 


-у/(5), є). - (0) дорівнює ге і максимумальна за умови Є| ,, --аа)/д, 


де «(а) - сила тяжіння на поверхні сферичної Землі. Припущення гармонічності сили тяжіння в 


перерахунку майже на висоту земного радіуса (6000 км) дає похибку - 245 Гал. 

Отже, для трансформації спостереженого (і нормального) поля сили тяжіння у зовнішньому просторі 
слід розв'язати |5| граничну задачу для рівняння сили тяжіння (2), конкретні оператори якого для різних ре- 
ференц-поверхонь |4| земного еліпсоїда обертання з малим стисненням не генерують істотно різних резу- 
льтатів поблизу Землі. Явні аналітичні вирази фундаментальних розв'язків рівняння (2) отримано в |ЄЇ. 

Похибки аномалій. Перерахунок ,аномалії" сили тяжіння до -- 2(5)- (5) (різниці нормаль- 
ного і спостереженого поля сили тяжіння) має свої нюанси |7). 

Через векторну природу до обчислене аномальне поле 


відрізнятиметься від істинної величини аномалії сили тя- 


ві«-ві- І НЯ -беоа де а - 


чено ее 


кут між нормальним ОЕ - У і виміряним 08 - є вектором 


жіння (рис. 1) 





сили тяжіння, СЕ - проекція вектора аномалії ВЕ на на- 
прям У, ЕР «| -|8| - звичайне значення аномалії, СР - 
його відхилення від проекції аномалії. Підставивши замість 


о це 2 
І рр, маємо ІвР)- А/Д 





п чан она 
уд 74 4 


Іза Ерб ЦИ -| 8-0) аномалія сили тяжіння може бути 





7 значною, оскільки ВЕ -2|мзіпац/2. і при й«їі5" ВЕ| «то 


МОЗ МЕОМЕРРИЧНЕ УМЕНЕ НА АНОМАЛІЙ мГал. На цю величину можуть відрізнятися звичайні анома- 


лії сили тяжіння Й -|8| від точних. 
Якщо б значення 7 -|8| давали високоточну аномалію сили тяжіння, її неможливо перера- 


хувати у зовнішньому просторі: рівняння сили тяжіння істотно залежить від швидкості зміни ку- 
тів між напрямом сили тяжіння (аномалії) і координатними осями. В локальних аномаліях цей 
напрям може сильно варіювати і будь-яке наближене припущення щодо нього істотно спотво- 


рить перераховані аномалії. Вираз 7-18, з точністю СР « Зо 8іп" с/2 співпадає з проекцією 
вектора аномалії на напрям вектора У (за Юнгом); при є -15" СР -2:107 мГал. 

У потенціалі И/ аномальних мас, які генерують аномалію сили тяжіння, його частинні похідні 
в будь-якому напрямі співпадають з проекціями аномалій сили тяжіння за цим напрямом, 
1-2 - дУУ / дх, ,де х, - спрямоване за напрямом вектора У . На площах у кількадесят км змі- 


ною напряму У нехтують, і аномалія 7 -|2| майже гармонічна. З точністю 0«єх«2004 мГал 


аномалії сили тяжіння на рівнині є значеннями вертикальної складової за напрямом модуля гра- 
дієнта. 

Перерахунок у нижній півпростір з відомою густиною с(Р) в рівнянні (2) нескладний. Але в 
практиці регіональної гравіметрії при перерахунку у бік мас нерідко на поверхні 5 відомі складо- 
ві ОМУ /дх, аномалії, а під поверхнею аномальні маси розташовані на невідомій віддалі від неї і 


слід гармонічно продовжити ОМ / дж, у бік аномальних мас до найближчих особливих точок поля. 


У такій постановці ця задача - істотно некоректна (має нескінченну множину розв'язків |41). 
Чисельне моделювання. Чисельне розв'язання нелінійної граничної задачі (3) у просторо- 
вому випадку здійснене за ітераційним процесом 


- ди з 
ли, 70, | У ай, 22 ем (5), (7) 
і-1 ОХ, 


"ГУ 





дей, ій, | - п-ії (п-1)-е наближення, аДХи ) - і-й спрямовуючий косинус градієнта (лп- 1)-го на- 


п-і1 
ближення й, |. З метою перевірки збіжності ітераційного процесу (10) розв'язано зовнішні гра- 


ничні задачі 


П 
си 05 вади, а й 
1 7, з, 


П Пп 


Ли - 0, (огад ДАНЕ ни и Б НИ СВО р 
г иї, (хоп кю РЯВУЛ МЕРЕ 







Зх чок 30 хіх -2Ах НОЮ 
ПРИ Б ВА Р ДАРА А В Є У, 

2 2 2 2 2 2 з, 2 
(ж - п) тот» (ж -н 1) тота (с -п) пре, (сп) ПрЕзРЕУ, 


де п-10, х, -5 км, У - горизонтальна площина. Точний розв'язок граничних задач (8) і (10) 


ДА щ П ЕЕ П а 
ож ма 2 М 2 2 2 2 2 2 | 
Ух ВХ їх (х -Л) вок (ж --Л) о ер. 


в'язку не має. На рис. 2 наведені точні значення градієнта й, і похідних задачі (8) у 30 точках 


хіх - их Ю (10) 


Ли-0), 





згадці, -Й 


дорівнює и - ) точного роз- 


рівномірної мережі поверхні 5 з кроком Дх, - ЇІ0 км на в області -25 км ЗХ, 5 25 км, - 25 км 
«хо2 15 км при грубих наближеннях напрямних косинусів - ац(и )- аг (и, )-0, сь(ц )- 1. У 
іншому варіанті розв'язку а (и ) «0.8, сь (и )-0, ази, ) «0,6. Як видно з фрагменту таблиці, 


градієнт и, і права частини граничної умови (8) співпали з точністю є 210" за 5 ітерацій у 


першому варіанті і за 6 - у другому. Модулі градієнтів точного і наближеного співпали з 
точністю 10". 

Кожну з граничних задач розв'язу- 
675 --142 Складові градієнта вали для різних перших наближень і 
ати, положень допоміжних точок |51, що ви- 
значають фундаментальні розв'язки. 
Рис. 2 відповідає випадку, коли допо- 


міжні точки взяті на площині х, «І. 


Число ітерацій у варіантах не переви- 
щувало 10. Хоча граничні значення 
точних і наближених градієнтів співпа- 
ли з точністю ай, вони самі спів- 
пали з меншою точністю через погану 
обумовленість взагалі задач з похід- 
ними в граничних умовах (для них, на 
відміну від задачі Діріхле, не діє прин- 
цип максимуму) і це слід врахувати в 
практичному моделюванні. 


Точка х -х, сх, -0, яка ство- 


рює поле у задачі (11), є однією з до- 
поміжних точок |51, що визначають 
фундаментальні розв'язки. Це пояснює 
бавовни заніоваюл пад зако ван зв бод ал 4 6 6 10 12 14 16 18 20 22 24 26 28 30 чималу точність визначення похідних 














0 4 8 42 16 20 24 28 ди/дх, за наближеним значенням ро- 
о 
ме и, ди, ди, ди, зв'язку. Якщо ця точка не входить до 
складу допоміжних точок, точність ви- 
дх дх, Ох, 


значення похідних гірша, особливо для 
234561 3,875 3,875 0,519 ди/дх 
275,455 5,281 4188 0,839 и АЙ | з 
313,623 9.262 3,087 1,236 дослідженні похибки розв'язку 
ззоз Цлбво | РТ альних розичанів рівномірно по 
313,58 9,262 3,085 1,234 р Ак енер 
5 | 275,596 5,283 41391 0,837 ташованих на площині Хх, - І, межею 
275,607 4,190 5,282 0,837 5 є площина х, -5. Ліва і п з 
рава час 
8 | 350166 5,281 5,281 р З , 8у пі до 
з рано Ціоаат | 8535 | за дн ен нн з Мр Ра 
сна наче он ла: інн похибка розв'язку и на 
440203 15,269 8,833 3,412 за У з 


для похідних -- 


абсолютна - 
350,065 8,582 8,984 1,717 похибки сягають 8096. Збільшення 


Фис 2. Розв'язання зовнішньої граничної задачі Алексідзе (8)і ЧИСЛа функцій до 49 знизило похибку 


, : і М я є і 
фрагмент відповідної таблиці даних, х 107. розв'язку й і похідних до 3-50, | 


ди/дх. - до 2096. Дальше нарощення 


до 64 збільшило похибку розв'язку і похідних до ГОМ, хоча права і ліва частини граничної 
умови (8) після 8 ітерацій співпали з точністю Є- -10? . Цей ефект пояснюють |4| немінімальністю 
систем фундаментальних розв'язків і малістю пре НИ ВО матриці коефіцієнтів для 
фундаментальних розв'язків. 

Коли напрям градієнта на всій границі 5 складає з внутрішньою нормаллю гострий кут, зов- 
нішня гранична задача (3) розв'язується із задовільною точністю і похибка зумовлена не пога- 
ною збіжністю ітерацій (14) (число ітерацій мале), а поганою обумовленістю граничної задачі з 
косою похідною. 


Гранична задача з косою похідною уздовж напряму |/ 


Аи «0, 9 (8 (11) 
ЗІ, у(5) 


за умови, що напрям зовнішньої нормалі М неортогональний напряму І (І в жодній точці гра- 





ниці С не співпадає з дотичною) 


ів (М,І)» 0 (12) 


х, є5 
не має відмінного від сталого розв'язку |8|. У граничних задачах (8-10) ця умова виконується на 
будь-якій скінченній частині площини У , оскільки поле створене точковим джерелом, а на всій 


нескінченній площині - ні, оскільки в нескінченно віддаленій точці площини 5 (м, Г) -й0. 
Ха. 0 


Для внутрішньої граничної задачі (11) і точкового джерела Р поза областю С умова (12) для 
достатньо гладкої границі не виконується, оскільки завжди існуватиме точка на 5, в якій 


радіус-вектор з початком в Р ортогональний до нормалі. Чисельні експерименти на сфері С 


виявляють погану обумовленість граничної задачі (11): граничні умови задовольняються доб- 
ре, а обчислені і точні значення градієнту и і його похідних - істотно різні. 


Розв'язок внутрішньої граничної задачі (8) для одиничного куба С (зовнішня нормаль має 6 
різних напрямів, початок координат (0,0,0) - на відстані 0,2 від нижньої кромки куба, 24 допо- 
міжні точки узяті на кубі зі стороною 1,4) за (10) при ац(и )- а (и )-0, аз(и, )- 1 такий: після 5- 


ї ітерації граничні значення задовольнялися ( Є «19о6), а наближення и і його похідні не мали 


спільного з точним розв'язком (/ -р/ Хі хі; рю і його похідними; нарощення числа ітерацій 


до 24 не збільшило точність наближень. У табл. 1 у 4-х точках границі, що є точками коллокації, 
дані точні і наближені значення градієнта і похідних. У решті точок значення функції через си- 
метрію співпадають з одним з наведених в табл. 2. Додали 25-й фундаментальний розв'язок, 
який співпадає з точним розв'язком задачі (8). Після о ї ітерації граничні умови та значення ро- 
зв'язку і похідних задовольнялися з точністю Є з- ен . Щоб розділити вплив нестійкості гранич- 


ної задачі (11) на сумарну похибку розв'язку задачі (8), вирішено задачу (11) для функції 
С -р/ х р; рю . Останні стовпці табл. 1 дають її розв'язок ий у тих же точках: похибка зу- 


мовлена нестійкістю задачі (14), а не повільною збіжністю ітерацій (7). 


Таблиця 1. Розв'язок внутрішньої граничної задачі (8) та що з косою похідною (11) для одиничного куба. 


2 КА ЛА А 
ленню винен вена в 
41989 Ї3 |2 6 | 30 | 18 | 0 | 8 | 13 | 7 | 2 | 8 


Збіжність ітерацій (5) для плоского випадку має особливості: розв'язок зовнішньої граничної 





задачі (4) з граничною умовою у(5)а а ху) є-1 і точним розв'язком ис пах ху для 
5 


о 0 
одиничного круга з центром в початку координат з напрямами а - со5(1-нФ), 2 -8іп(г--ф) п 
- центральний кут, Ф - збурення початкового напряму градієнта шуканої функції) має похибки, 
що залежать від положення точок ітерації. У табл. 2 дано кількість ітерацій М до умови їх зупин- 





ки | - |вгайн| -1 «107, максимальні похибки розв'язків після кожної ітерації у точках колло- 


вин " -4 " " вини " у (я з4 
кації ДЛЯ ПОХІДНИХ - Є «10 , а посередині між точками коллокації в кінцевому розв язку Є» «10 ; 


и ди") ди 
де Є 2 Ро вна . 
Х Ж. Х; 


і і 


Таблиця 2. Параметри розв'язку плоскої зовнішньої граничної задачі (4) для одиничного круга. 





Розв'язок внутрішньої граничної задачі за сприятливих умов визначається з точністю до ста- 
лої (тому не порівнюємо точний і наближений розв'язки) і знаку (при 2-7 1-ша ітерація задо- 
вольняє граничній задачі). 


Для тієї ж задачі (4) при радіусі кола г -2 узяті числа 7 -д -2/2 (напрям градієнта 


співпав з прямою Х, - Х,). Цей напрям у точках У2 2 і -У2 92 кола утворює з нор- 
2 2 7, ПУ 

маллю кут 907, тому не виконується умова (12). Точний розв'язок цієї задачі и - 0.693, а значення 

граничної умови у/(5) - 0,5. Рис. З ілюструє розв'язок задачі (4) (9 ітерацій, Х, Х» - координати 

граничних точок коллокації). Градієнти наближеного і точного розв'язку близькі, самі розв'язки і 

їхні похідні істотно різні, оскільки наближення 7 -д) -2/2 не забезпечує умови (12). 






-1.5 м 2.5 о 05 у -1 з -0, 05 


Рис. 3. Розв'язок зовнішньої граничної задачі (4) для кругаз Г -- Й при 7 -д) -/2/2 


Отже, для збіжності ітерацій (5) необхідно не лише, щоб точний розв'язок задовольняв умові 


(12), а і вибір початкових наближень, для яких в ітераціях (5) не виникали б коефіцієнти СІ; , ЯКІ 


не забезпечують умови (12). 
Розв'язання внутрішньої граничної задачі для одиничного круга з центром в початку коорди- 


нат і граничною умовою у(8)-/ (о. -17 (Се, -1) за ітераціями (5) здійснено за напрямом 
градієнта 7 - 0.8315, 5 - 0.5556, що, як і будь-який постійний напрям, не забезпечує виконан- 
ня умови (12). Рис. З ілюструє наближений розв'язок задачі після 8 ітерацій. Точний розв'язок 


2 2 | | 
ис іп (ж -1) --(х, -1) , але комплексно-спряжена ( Фаи-ніу - голоморфна функція) з ним 


функція у завдяки відношенню Коші-Рімана теж задовольняє цій задачі. За цих 7 ітерації (5) 


збігаються не до точного розв'язку и, а комплексно-спряженої функції у, що особливо видно 
на прикладі горизонтальних похідних (рис. 4). 





Рис. 4. Розв'язоквнутрішньої граничної задачі (4) для одиничного круга 


Моделювання аномалій. На практиці відшукують не потенціал у -и-у сили тяжіння, а по- 
тенціал и аномалії сили тяжіння, де потенціал нормальної сили тяжіння у вважають відомим. 
За лінійних наближень граничних задач граві- і магнітометрії цей перехід до аномалії не змінює 
не лише рівняння, а і вигляд оператора граничних умов, необхідно лише за граничну функцію 


у (5) узяти аномалію сили тяжіння. 
У нелінійних граничних задачах перехід до аномалії змінює ліву частину граничних умов (3, 
4), так як операція обчислення модуля градієнта не дистрибутивна: гади |вгайи| - |егад у. Як- 


що гранична функція у (5) - модуль градієнта потенціалу у/ сили тяжіння, для потенціалу ано- 
малії и одержуємо нелінійну граничну задачу: 


2 
3 д В 
Р и бу ди й 
Ди-0, УЖ ра -2У; рн (5)-|огаду 

ОХ Р ОС 
граничні значення якої за нульового потенціалу нормальної сили тяжіння співпадають з гранич- 
ними умовами (3). Для розв'язання цієї задачі задіяно один з варіантів таких послідовних набли- 
жень (умови на 00 опущені ): 


 ирву (13) 


ор 
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1 і 5 
З'ясування чисельної швидкості збіжності ітерацій (14) має практичний інтерес, бо на кожній 
ітерації розв'язується гранична задача і слабкозбіжні процеси істотно погіршують практичну екс- 
плуатацію алгоритмів. 


Збіжність (14) обчислено методом |4|) на границі 5 у вигляді площини Х3 - 5 на моделі 


РОН ХВ У РН З 


2 2 2 
уд ТТ дача и) 
на 36 вузлах коллокації в точках (Хі - -25410К, Х» « -25410), хз - 5, К, | - 1,2, 3, 4, 5) границі 5 
для 36 фундаментальних розв'язків, полюси яких є під вузлами коллокації на площині пана 


Ітерації тривали до виконання умови тах ||егаф(у ни, | -м(5 | «є: При єЄ:107 1-й ітераційний 
5 
процес завершився за 3 ітерації, другий - за 36, а останній - за 103 ітерації, а при є210"- за 
4, 79 | 147 відповідно. 
Для першої з ітерацій (14) розв'язана гранична задача з похідною уздовж Х. на границі. До- 
поміжні точки для фундаментальних розв'язків узяті на площині х, - 1. Наближення розв'язку і 


його похідних добре співпадають між собою (при є «107 з точністю 1072, при є -10" з точністю 
107), але далекі від точних. У табл. 3 дані їх середньоквадратичні значення у вузлах коллокації 
(1-й стовпець) і у середині між ними (2-й стовпець), а в табл. 4 - відповідні похибки для різних 
є-10"Ч(ка21,2..,5Уї М. 


Таблиця 3. Наближені розв'язки нелінійної граничної задачі (14) на площині. 


бі 
0.598 0.029 0.021 0.029 0.021 0.026 0.067 
0.611 0.032 0.028 0.032 0.028 0.038 0.053 
0.614 0.030 0.032 0.030 0.032 0.047 0.044 


бамідк, 


0.025 
0.018 
0.018 
0.026 
0.01 
0.026 
0.003 | 0.003 
0.002 | 0.002 

За табл. 4 виявлено, що збільшення точності апроксимації граничних умов у вузлових точ- 
ках коллокації без збільшення кількості вузлових точок (їі, отже, збільшення числа функцій, що 
беруть участь у розкладанні) не впливає на похибку (починаючи з К - 2). Ці обставини виявляють 
нелінійну поведінку похибок чисельного розв'язання граничних задач з похідними і спонукають досліджува- 
ти похибки розв'язків кожної конкретної задачі. 

Висновки. Трансформацію значень поля сили тяжіння у зовнішньому просторі можна здійснити через 
розв'язання граничної задачі (3) для рівняння сили тяжіння (2), конкретний вигляд якого залежить від обра- 
ної моделі Землі. У плоскому випадку ця операція зведена до розв'язання граничної задачі (4) з ітерацій 
(5). Точність її розв'язання істотно залежить від міри обумовленості задачі, наближень напряму шуканого 
градієнта, параметрів чисельного методу |5| (числа і положення фундаментальних розв'язків), та напряму 
зовнішньої нормалі, уздовж якої обчислюють похідні в граничних умовах (виконання умови (12). Попри це, 
її розв'язання на простих моделях за певних обмежень - успішне. 

Основними джерелами неоднозначності розв'язків задачі є скінченна вертикальна роздільна 
здатність розв'язку, неадекватний вибір моделі задачі чи початкових наближень розв'язку. 


0.056 
0.057 
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Точнье трансформации гравитационного поля в региональном масштабе обеспечивает реше- 
ние нелинейной граничной задачи Алексидзе для потенциала силь тяжести. Ее замена на соот- 
ветствующую задачу для уравнения Лапласа генерирует погрешности, оценки которьктх приведе- 
нь в статье. Численное моделирование последней задачи показало, что точность ее решения су- 
щественно зависит от мерьг обусловленности задачи, приближений направления искомого градивента, па- 
раметров численного методи и направления внешней нормали, вдоль которой вьтшчисляют производньке. 


